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Solvent  effects  on  electron  transfer  reactions  have  been 
the  subject  of  considerable  interest  in  recent  years.  For  a 
simple  heterogeneous  process  such  as 

A  +  e  B  (1) 

where  both  A  and  B  are  present  in  an  electrolyte  solution,  one 
may  recognize  four  effects  of  the  solvent  on  the  kinetic 
parameters.  Firstly,  the  standard  potential  for  the  reaction 
depends  on  the  free  energy  of  solvation  of  reactant  and  product, 
and  changes  with  the  charges  on  these  species  and  the  extent  to 
which  a  given  polar  solvent  can  solvate  the  reactant  and  product 
species  [1].  The  second  way  in  which  the  solvent  plays  a  role  is 
through  the  double  layer  effect  according  to  which  the 
concentration  of  the  reactant  at  the  reaction  site  and  the 
potential  drop  experienced  by  the  electron  during  transfer  are  a 
function  of  the  local  electrostatic  field,  which,  in  turn, 
depends  on  the  charge  density  on  the  electrode  and  the  dielectric 
properties  of  the  electrolyte  solution  in  the  double  layer  [2]. 
Thirdly,  the  solvent  affects  the  free  energy  of  activation  for 
reaction  (1),  especially  through  the  dependence  of  the  outer 
sphere  contribution  to  solvent  reorganization  on  the  dielectric 
properties  of  the  surrounding  medium  [3].  This  aspect  has  been 
of  particular  interest  in  recent  years,  many  workers  having 
examined  the  validity  of  the  elegant  though  perhaps  simplistic 
Marcus  dielectric  continuum  treatment  [4-8].  Finally,  the 


dynamic  properties  of  the  solvent  determine  the  frequency  which 
the  activated  complex  is  formed,  and  further,  the  frequency  with 
which  such  a  complex  successfully  decays  to  the  product 
configuration.  This  aspect  is  expressed  in  the  pre-exponontial 
part  of  the  standard  rate  constant,  and  was  considered  in  detail 
in  two  recent  seminal  papers  by  Calef  and  Wolynes  [9,10].  since 
the  appearance  of  their  work,  a  number  of  experimental  studies  of 
solvent  effects  on  heterogeneous  electron  transfer  reactions 
involving  organic  and  organometallic  molecules  have  been  reported 
in  the  literature  [11-15].  The  systems  chosen  were  ones  in  which 
the  component  of  the  free  energy  of  activation  due  to  internal 
reorganization  was  small  so  that  the  effect  of  the  solvent  on  the 
observed  standard  rate  constant  was  easily  seen.  Particularly 
impressive  is  the  work  of  Opallo  [14]  who  studied  the 
electrooxidation  of  1, 4-diaminobenzene  in  twelve  different 
solvents  and  obtained  kinetic  parameters  in  a  range  where  good 
precision  is  possible.  Opallo  also  demonstrated  in  this  work, 
and  in  an  earlier  study  with  Kapturkiewicz  [13]  that  solvent 
effects  are  clearly  different  for  protic  and  aprotic  solvents. 
Nonetheless,  a  clear  correlation  between  solvent  dynamical 
properties  and  the  standard  rate  constant  was  demonstrated. 
Similar  results  were  obtained  by  Weaver  and  coworkers  [12,15], 
who  examined  the  reduction  kinetics  of  various  metallocene 
complexes  and  iron  cyclooctatetraene  tricarbonyl  in  a  series  of 
aprotic  and  protic  solvents. 

At  first  inspection,  the  results  of  both  studies 
mentioned  above  show  no  indication  of  the  connection  between  the 
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rate  constant  and  the  solvent  reorganlzational  energy  estimated 
from  the  Marcus  treatment.  This  result  is  surprising  given  the 
excellent  agreement  between  absorption  maxima  for  intramolecular 
charge  transfer  in  mixed  valence  complexes  and  the  Marcus 
predicted  outer  sphere  reorganization  energies  [4].  In  the 
present  communication  we  consider  in  detail  the  general 
expression  for  the  rate  constant  proposed  by  Calef  and  Wolynes 
[9,10],  and  present  an  analysis  that  demonstrates  the  validity  of 
their  form  of  the  pre-expontial  factor  as  well  as  the  Marcus 
dielectric  continuum  treatment  for  the  activation  free  energy. 
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^  Analysis 


According  to  the  encounter  pre-equilibrium  model  [9,10], 
the  standard  rate  constant  corrected  for  double  layer  effects  is 
given  by 


k,c  =  i:KpKnexp(-AG*/RT) 


(2) 


where  k  is  the  electronic  transmission  coefficient,  Kp,  the 
equilibrium  constant  for  precursor  complex  formation,  the 
nuclear  frequency  factor,  and  aG*,  the  reorganizational  free 
energy  of  activation.  The  latter  quantity  is  made  up  of  an  inner 
sphere  contribution  aG*  j  g ,  which  may  be  calculated  from  the 
vibrational  force  constants  for  the  reactant,  and  the  coordinate 
changes  accompanying  electron  transfer  [16,17],  and  an  outer 
sphere  contribution  ^G*g8  given  by  [3] 


aG* 


0  s 


Noe* 

32nfo 


(3) 


where  Ng  is  Avogadro's  number,  e,  the  charge  on  the  electron,  a, 
the  radius  of  the  reactant  represented  as  a  sphere,  R,  the 
distance  of  the  charge  center  of  the  reactant  from  its  image  in 
the  electrode,  fgp,  the  optical  dielectric  constant,  the 
static  dielectric  constant,  and  the  permittivity  of  free 
space.  This  estimate  of  aG*o,  is  based  on  the  Born  model,  and 
inherits  some  key  simplistic  assumptions,  namely,  that  the 
solvent  is  a  structureless  continuum  and  that  there  is  no  spatial 


-5- 


dispersion  of  the  dielectric  constant  at  distances  beyond  the 
reactant  radius. 

For  systems  in  which  aG*  o  s  i  s »  nuclear 

frequency  factor  is  described  by  the  overdamped  solvent 
relaxation  model  and  is  given  by  [9,12] 


n 


1 


.  4itRT  . 


1/2 


(4) 


where  Ti  is  the  longitudinal  solvent  relaxation  time.  This 
quantity  can  be  determined  from  the  dielectric  properties  of  the 
solvent  as  follows: 


Tl  -  -  (5) 

*  s 

Tq,  is  the  usual  Debye  relaxation  time  and  the  dielectric 
constant  of  the  solvent  at  infrared  frequencies.  Values  of  Tl 
have  been  estimated  for  a  large  number  of  aprotic  solvents 
commonly  used  in  non-aqueous  electrochemistry  [13,15]  and  are 
summarized  in  Table  1.  It  is  clear  that  there  is  considerable 
variation  in  t^,  namely,  from  0.2  ps  in  acetonitrile  (AN)  to  8.9 
ps  in  hexamethylphosphoramide  (HMPA) .  It  is  important  to  point 
out  that  eq.  (4)  is  strictly  applicable  only  to  solvents  for 
which  Tl  is  greater  than  l  ps  [9].  Data  discussed  in  this  paper 
certainly  include  solvents  that  do  not  meet  this  criterion,  a 
point  which  is  examined  further  below. 
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On  the  basis  of  the  models  described  here,  it  follows 
that  the  corrected  standard  rate  constant  kg^  depends  on  the 
nature  of  the  solvent  in  two  ways,  namely,  through  variation  of 
the  nuclear  frequency  factor,  v^,  and  variation  in  the  outer 
sphere  contribution  to  the  free  energy  of  activation,  aG*os. 

When  one  examines  the  data  published  to  date  [11-15],  it  is  clear 
that  variation  due  to  the  nuclear  frequency  factor  predominates. 
Thus,  Opallo  [14]  obtained  an  excellent  linear  correlation 
between  the  logarithm  of  the  standard  rate  constant  and  the 
logarithm  of  the  longitudinal  solvent  relaxation  time  Ti  using 
data  for  the  oxidation  of  1,4-diaminobenzene  in  eight  different 
aprotic  solvents,  in  spite  of  the  fact  that  is  changing  at 

the  same  time.  This  result  is  not  surprising  when  one  considers 
the  fact  that  changes  by  a  factor  of  40  for  the  solvents 
considered  by  Opallo  [14]  whereas  the  solvent 

dependent  part  of  changes  only  by  a  factor  of  1.4.  In 

addition,  there  is  no  correlation  between  and  (e’^op  "  ^'^s) 
(see  Table  1) .  In  analyzing  their  data,  previous  authors  have 
attempted  to  confirm  the  validity  of  the  encounter  pre¬ 
equilibrium  model  [9,10]  by  estimating  values  of  kgj  on  the  basis 
of  the  above  equations  and  comparing  the  estimates  with 
experimental  results.  This  procedure  is  rather  unsatisfactory 
because  it  requires  that  one  assume  values  for  the  parameters  a 
and  R  in  order  to  estimate  aG^o,,  and  a  value  for  rKp  in  order  to 
estimate  the  rate  constant,  kg^.  Weaver  [12,3.5]  compared 
calculated  and  observed  values  of  kg^  by  plotting  the  logarithm 
of  one  quantity  against  the  logarithm  of  the  other.  In  such  a 


comparison,  the  choice  of  cKp  affects  only  the  intercept  of  the 
plot  whereas  the  value  chosen  for  the  distance  parameter  aR/(R-a) 
affects  the  slope  and  the  correlation  coefficient  in  a  least 
squares  analysis. 

It  seemed  to  us  that  these  data  provided  an  opportunity 
to  test  the  validity  of  the  Marcus  formulation  [3]  for  -^0*^8  in 
so  far  as  its  solvent  dependence  is  concerned.  Combining  eqs. 
(2)-(4),  and  defining 


op 


(6) 


and 


Noe‘ 


32it«8RT 


'1  V 

.a  R/ 


(7) 


one  obtains  the  expression 


^8  C 


tKp  fg^^ 


1/2 


\4Tr/ 


exp(-AG*i ,/RT)oxp(-g5) 


(8) 


Rearranging,  one  may  write 


In 


fkscU) 


=  In 


tKp’ 


4ivj 


RT 


-  g5 


(9) 


It  follows  that  a  plot  of  ln(k,jTL/8^'^*)  against  8  should  be 
linear  with  a  slope  of  -g.  Furthermore,  if  aG* j ,  is  known,  one 
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may  es^.mate  the  product  cKp  from  the  Intercept.  Thus,  the 
present  analysis  assumes  that  the  dependence  of  the 
preexponential  factor  on  the  solvent  parameters  Ti  and  s  is 
correct  as  described  by  the  encounter  pre-equilibrium  model 
[9,10],  and  examines  the  dependence  of  the  kinetic  data  corrected 
for  part  of  the  solvent  effect  on  the  dielectric  parameter  8  in 
order  to  assess  the  validity  of  the  Marcus  expression  for  the 
outer  sphere  contribution  to  the  reorganizations!  free  energy. 
Details  of  our  analysis  of  data  from  the  literature  [11-15]  are 
now  discussed  in  detail. 


*This  equation  was  used  previously  by  Harrer  et  al.  [23]  to 
demonstrate  solvent  effects  on  the  rate  constant  for  homogeneous 
electron  transfer  between  tetracyanoquinodimethane  and  its  anion 
radical.  However,  these  authors  did  not  analyze  the  parameters 
obtained  from  the  linear  fit. 
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Results  and  Discussion 


The  above  analysis  was  applied  to  four  sets  of  kinetic 
data  for  systems  In  which  the  kinetic  parameters  had  been 
obtained  In  at  least  four  aprotlc  solvents  using  the  same  base 
electrolyte  [13-15].  In  all  cases,  a  good  linear  correlation  was 
obtained  between  the  solvent  corrected  kinetic  parameter 
ln(kg  5  Tl/8  ^ ^ )  and  the  permittivity  parameter  5  with  a  negative 
slope.  While  previous  work  has  demonstrated  that  the  solvent 
exerts  its  influence  on  electron  transfer  primarily  through  the 
nuclear  frequency  factor  the  present  analysis  makes  it  clear 
for  the  first  time  that  the  Marcus- expression  for  is  also 

valid  for  heterogeneous  electron  transfer. 

Values  of  the  slope  and  intercept  for  plots  of  published 
data  [11-15]  are  summarized  in  Table  2;  in  addition,  the  plots 
obtained  for  the  1, 4-diaminobenzene  [14]  and  cobaltacenium 
systems  [15]  are  shown  in  Figures  1  and  2,  respectively.  In  all 
cases,  the  magnitude  of  the  slope  yields  very  reasonable  values 
for  the  reactant  size-distance  parameter,  aR/(R-a).  However,  the 
values  obtained  here  are  quite  different  from  the  values  assumed 
by  Opallo  [14]  and  Weaver  et  al.  [15].  For  instance,  in  his 
analysis,  Opallo  took  the  reactant  radius  a  for  1,4- 
diaminobenzene  to  be  0.34  nm,  as  determined  from  diffusion 
coefficient  data  using  the  Stokes-Einstein  equation.  Assuming 
the  imaging  distance  R  to  be  infinity,  the  value  of  aR/(R-a)  is 
0.34  nm  which  compares  with  the  value,  0.62  nm  obtained  by  the 
present  analysis.  Similarly,  Weaver  et  al.  [15]  assumed  a  -  0.38 
nm  and  R  ■  «  in  analyzing  data  for  the  cobaltacenium  system. 
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According  to  the  present  analysis  aR/CR-a)  Is  equal  to  0.85  nm,  a 
considerably  higher  value.  These  results  are  significant  In  two 
respects.  Firstly,  In  selecting  values  for  the  parameter  aR/(R- 
a) ,  previous  authors  have  essentially  doubled  the  magnitude  of 
^G*os*  Secondly,  If  one  calculates  R  on  the  basis  of  their 
assumed  molecular  radii  a,  which  are  certainly  reasoneUsle,  It  Is 
apparent  that  R  la  finite  and  the  order  of  0.7  nm.  Indeed,  the 
assumption  that  R  >  «  Is  a  precarious  one  theoretically,  since 
the  electronic  transmission  coefficient  approaches  zero  for  large 
values  of  R,  and  the  reaction  can  no  longer  be  considered 
adiabatic.  In  any  case,  some  disagreement  between  the  true  and 
calculated  values  of  aR/(R-a}  Is  not  unexpected  considering  the 
fact  that  the  Marcus  model  assumes  simple  spherical  Ions  are 
Involved  In  the  reaction.  Thus,  an  Improved  estimate  of 
can  be  obtained  If  a  species  such  as  the  cation  radical  formed 
from  1, 4-dlamlnobenzene  Is  represented  as  a  collection  of  spheres 
rather  than  a  single  sphere  [18]. 

For  both  the  1, 4-dlamlnobenzene  and  cobaltacenlum 
systems,  estimates  of  have  been  made  so  that  values  of  cKp 

could  be  extracted  from  the  Intercepts  of  the  plots  shown  In 
Figs.  1  and  2.  For  1. 4-dlamlnobenzene,  Grampp  and  Jaenlcke  [19] 
estimated  aG*|,  to  be  0.9  kJ  mol~^.  Accordingly,  from  the  value 
of  the  Intercept,  the  estimate  of  cKp  Is  0.48  ±  0.48  pm.  While 
the  error  In  the  estimate  Is  large,  cKp  Is  still  two  orders  of 
magnitude  smaller  than  the  value  of  60  pm  proposed  by  Hupp  and 
Weaver  [20]  and  assumed  by  Opallo  [14].  For  the  cobaltacenlum 
case,  when  ^^G*  ^ ,  Is  estimated  to  be  1.0  kJ  mol~^  [15],  cKp  Is 
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2.7  ±  2.8  pa.  While  larger  than  that  for  the  1, 4-diaminobenzene 
case,  this  result  also  suggests  that  previous  estimates  of  tKp 
are  too  high. 

In  assessing  the  above  results  one  should  keep  in  mind 
the  fact  that  the  longitudinal  relaxation  time  Ti,  which  has  the 
major  influence  on  the  observed  rate  constants,  varies  with  the 
nature  and  concentration  of  the  supporting  electrolyte  [21-22]. 
The  values  cited  previously  [13-15]  and  summarized  in  Table  1  are 
those  for  the  pure  solvent.  On  the  basis  of  the  kinetic  data 
obtained  by  Kapturklewicz  and  Opallo  [11]  for  the 
electroreduction  of  nitromesitylene,  it  is  clear  that  the 
standard  rate  constant  after  double  layer  correction  varies 
significantly  with  the  nature  of  the  cation  of  the  electrolyte. 
Thus,  the  value  of  kgj  obtained  in  tetrabutylammonium  perchlorate 
was  always  smaller  than  that  found  in  tetraethyl ammonium 
perchlorate.  This  observation  can  be  attributed  to  a 
corresponding  decrease  in  the  viscosity  of  the  electrolyte 
solution,  and  thus,  in  the  value  of  t^.  Considering  the  fact 
that  few  dielectric  relaxation  data  are  available,  it  is  clear 
that  one  should  only  attempt  the  type  of  analysis  presented  here 
for  systems  in  which  the  concentration  and  nature  of  the  base 
electrolyte  is  kept  constant  when  the  solvent  is  changed.  Then, 
the  assumption  made  is  that  the  actual  values  of  parallel 
those  observed  in  the  pure  solvent,  and  that  the  difference 
between  these  quantities  affects  only  the  values  of  the  intercept 
of  plots  such  as  those  presented  here.  It  should  also  be 
remembered  that  the  values  of  Ti  recorded  for  some  solvents  are 
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rather  Inpreclse.  Although  reasonably  precise  values  of  the 
Debye  relaxation  time  Tq  are  available  for  most  solvents,  the 
values  of  are  sometimes  uncertain  due  to  the  fact  that  the 
frequency  range  for  the  dielectric  measurements  was  not 
sufficiently  high.  The  availability  of  more  data  for  dielectric 
relaxation  in  electrolyte  solutions  would  do  much  to  improve  the 
situation. 

It  should  also  be  pointed  out  that  in  both  the  present 
analysis  and  previous  work  [11-15],  it  has  been  assumed  that  the 
image  distance  R  is  independent  of  the  solvent.  If  the  reaction 
site  is  on  the  outer  Helmholtz  plane,  this  assumption  is  strictly 
not  valid  since  the  thickness  of  the  inner  layer  depends  on  the 
size  of  the  solvent  molecule.  This  is  obviously  connected  to  the 
other  assumption  implicit  in  this  analysis,  namely,  that  cKp  is 
also  a  solvent  independent  parameter.  The  electronic 
transmission  coefficient  c  is  certainly  distance  sensitive,  and 
hence,  solvent  sensitive  as  well  if  the  reaction  site  is  on  the 
outer  Helmholtz  plane  or  fuxrther  out  in  the  diffuse  layer,  with 
the  present  data,  however,  we  can  find  no  connection  between  the 
molecular  volumes  for  a  given  solvent  and  deviation  from 
linearity. 

Finally,  as  a  comment  on  the  general  validity  of  the 
overdamped  solvent  relaxation  expression  for  we  note  the  fact 
that  some  of  the  solvents  in  which  longitudinal  relaxation  is 
rapid,  namely,  acetonitrile,  dichloromethane,  and  tetrahydro- 
furan,  to  which  eg.  (4)  is  not  strictly  applicable,  seeing  to  fall 
nicely  on  line  with  the  present  analysis.  Perhaps  further  work 
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will  determine  whether  this  is  due  to  the  effect  of  the 
supporting  electrolyte  on  Ti,  or  to  a  fortuitous  shortcoming  of 
the  theory. 


We  are  presently  considering  solvent  effects  on 
homogeneous  outer  sphere  electron  transfer  reactions.  Meanwhile, 
it  is  hoped  that  the  present  work  will  help  stimulate  future  work 
and  solidify  experimental  strategies  for  the  study  of  solvent 
effects  in  heterogeneous  electron  transfer  reactions. 
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Values  of  tj  and  are  from  the  compilation  of  Marcus  [1]-  Values  of  e„  and  Tq  were  taken  from  the  papers  of  Opallo  [14]  and 
Weaver  et  al.  [12,  15]  who  cite  the  original  literature. 


Table  II.  Summary  of  Parameters  Found  in  Fitting  Equation  (9)  to  Kinetic  Data  for  Heterogeneous  Electron 
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Legends  for  Figures 


Figure  1. 


Figure  2. 


Plot  of  the  solvent  corrected  kinetic  parameter 
In  against  the  permittivity  para¬ 

meter  l/«op  "  V^s  using  data  for  the 
electrooxidation  of  1,4-diamlnobenzene  at 
platlnxim  in  seven  aprotic  solvents  [14]. 

Plot  of  the  solvent  corrected  kinetic  parameter 
In  against  the  permittivity 

parameter  l/fop  “  V^s  using  data  for  the 
electroreduction  of  the  cobaltacenium  cation 
in  eight  aprotic  solvents  [15]. 
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